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Abstract 

In this paper, we prove a result of equivalence in law between a diffusion conditioned with respect 
to partial observations and an auxiliary process. By partial observations we mean coordinates (or 
linear transformation) of the process at a finite collection of deterministic times. Apart from the 
theoritical interest, this result allows to simulate the conditioned diffusion through Monte Carlo's 
method, using the fact that the auxiliary process is easy to simulate. 
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1 Introduction 

We are interessed in multidimensional diffusions solutions of stochastic differential equations (SDE's) 
generated by a Brownian motion. For a n-dimensional diffusion solution on [0, T] of the following 

Axt ^bt{xt)dt + at{xt)dwt., xq = u (1) 

where w is a n-dimensional Brownian motion, it is known (see e.g. [7]) that its conditional law 
^{x\xt — v) is given by the law of a bridge process (as extension of Brownian bridge) y solution of 

dyt = bt{yt)dt + at{yt)dwt + at{yt)crt{yt)*^z'^ogpt^T{z,v)\^^_y^dt, yo ^ u 

where ?i is a Brownian motion and Ps^t{z, ■) is the density of Xt knowing Xg = z. But in most cases this 
density is not explicitely known so that we are not able to simulate it easily. For practical purposes, e.g. 
parameter estimation of diffusion processes, simulation of paths corresponding to the conditional law is 
needed. 

In their paper [4], B.Delyon and Y.Hu studied the following equation on [0,r] 

dyt ^ btiyt)dt - j^—^dt + atiyt)dwt, yo = u. (2) 

where is a n-dimensional Brownian motions. Under adequate assumptions, the process y is unique on 
[0,T], \imt-^Tyt — V, a.s. and for all positive function / in C([0,T],R") we have 

E[f{x)\xT^v]^CE[f{y)R{y)] 

where i? is a functional of whole path y on [0, T]. The quantity R{y) is computable knowing parameters 
b, a, T and v. The constant C is unknown, but in practice the conditional law is estimated through 

EM 

where each is an independant sample of ([1]). In this case, we call the process y a bridge even if y does 
not have the right targeted law. If 6 = and a — In (identity n-dimensional matrix) , the process a; is a 
n-dimensional Brownian motion and process y is a 7i-dimensional Brownian bridge so that C = R = 1. 
This theorem applies in the case of more than one observation. The Markov property indeed implies 
that the conditional law is the tensor product of each bridge. 
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The aim of this paper is to extend this result to solve this problem with only partial observations. 
The previous remark does not apply; indeed we have to treat simultaneously all conditionings . To give 
an idea, let w = ( ^2 ) be a 2-dimensional Brownian motion. The law of w conditioned on Wg — u and 
w'^ = V with 5 < T is given by that one of y solution of 

Szl '^"^ dt, yo = u 

each coordinate is a Brownian bridge. 

Let us define our observations. At each deterministic positive observation time of the sequence 
< Ti < • • • < Tfc < • • • < Tjv = T, we get a partial information given by a linear transformation of xt^. , 
LkXTf., where is a deterministic matrix in M^^, „(R) whose vrik rows form an orthonormal family. 
So that our aim is to be able to describe the conditional law J^{x\{LkXTk = Vk)i<k<N) where Vk is an 
arbitrary deterministic rrifc-dimensional vector. 

We define process y to be the solution of 

\ yo = u 

where for all time t, all vector z and for all 1 < A: < TV, the matrix Pt{z) is an oblique projection and 
Uk is any vector satisfying LkUk = Wfc. The correction term operates only on the interval (Tk — ek,Tk) 
where Tk ~ Sk < Tk for technical reasons. We will show that with a good choice for those projections 
(see Equation ([5])) we have the following equivalence in law 



^{x\{LkXT^ = Vk)i<k<N) ^ ^{y) 



with an explicit density (Theorem [T] below). 

In this paper a first part is devoted to the study of general bridges which will provide us the good 
candidate whose law is absolutely continuous with respect to targeted one. The second one provides the 
main result. Some properties and proofs are postponed in the appendix to ease the reading. 

Notations For the sake of readibility, we choose not to specify arguments when not necessary. For 
example ((!]) becomes 

dxt — bt.dt + atdwt 

For all z, the matrix at{z) is defined by 

at{z) = at{z)at{z)* 
we suppose that there exists a positive number p such that for all {t, z) 

p^^In < at{z) < pin 

in the sense of symmetric matrices, where /„ is the n- dimensional identity matrix. The function is 
defined by 

: [0,r]xR" ^ A/„(R") 

We define the infimum of all the Sk 

£0 = minkiek} 
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2 Bridges and bridges approximations 

Bridges 

We recall that a bridge is defined as a solution of ([3]) 

VO = Mo 

We assume that the deterministic parameters b and a are C^'^ functions (bounded with bounded deriva- 
tives). We assume that 

is a C^'^ function and that for any z 

LkPt\z)=Lk and ker(Lfc) - ker(P,'=(z)) (4) 
First of all, a lemma to describe the behaviour of process y 

Lemma 1. The SDE admits a unique solution on [0,T] in the absolute convergence's sense meaning 
that 

\\PHyt){yt-uk)\\ ^,^ ^ ^ 

at < +00 

For all k we have Lkyr^ = LkUk almost surely (a.s.). Moreover for Tk — Sk < t < Tk, \\Lk{yt — Uk)\\ < 
Ck{oj){Tk — t) loglog[(rfc — t)^^ + e] a.s., where Ck is a positive random variable. 

Proof. Let us remark that for times in [Tk-i,Tk] (with Tq — 0) the SDE becomes 

dyt = btdt - P^yi-J^l^^^_^^^^^){t)dt + aiyt)dm 

So that we may reduce the proof to the study of ^ with only one observation time, but we here have 
to consider random initial conditions. If unicity holds it will lead to the result by concatenation. The 
proof in the case iV = 1 is given in the appendix with Lemma ID □ 

Bridges approximations 

We now introduce approximations that will be useful in the proof of the main result in next section. Let 
< e < So, we set 

dyl = bl{yl)dt + at{yt)dwt - ^ P^\yl)ylz^i^^^^^^^^^_^^(t)dt, yl = (5) 

k 

The only difference with the Bridge Equation Q is that each correction term is stopped from a distance 
e from the observation time. 

Lemma 2. There exists a constant < k < 1 such that 

sup EM-ytf]<Ce^ 
te[o,T] 

for all < e < {eq A 1) where C is a positive constant. The numbers C and k depend on T , N , the (sk)k, 
the {Ak)k, and the bounds for b and a. 

Proof. Given in the appendix, the proof uses classical techniques and auxiliary processes each defined 
on[Tfe_i,Tfe]. □ 
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3 Result in the case of partial observation 

Case where b is bounded 

We aim to obtain a Delyon&Hu-type theorem that gives absolute continuity of process x sohition of ([1]) 
conditioned on observations (LkXr^ = Vk)i<k<N with respect to a bridge process y solution of ([3]). 
We now consider a peculiar projection P, for all k and z 

Ptiz) ^ at{z)Ll{Lkat{z)Lir^Lk (6) 



We set 
and also 

Let us remark that 



A^iz) = {Lkat{z)Ll)- 



Pt{zf = at{zyLlA^^{z) and r^k{z) = ^det(A^(z)) 



/3f(z)*/3f(z) = Af(z) and L,a,(z)/3f (z) - (7) 
where Im^ is the identity m^— dimensional matrix. Here are both systems we now consider 

dxt = bt{xt)dt + (7tixt)dwt, xq ^ u 

dyt = bt{yt)dt + at{yt)dwt-^<7t{yt)l3t{yt)^^;^—j—l{n-e,,n)dt, yo = u (8) 



Theorem 1. Suppose b, a and a ^ to be Cl'^ -functions. Then for any bounded continuous function f 



The result is the following 

E[/(a;)|(Li;XTfc ^ Vk)i<k<N] 



= CE /(y)_[J_77fe(yTjexp| ^-^ — +/ — — 

{Luys - Vk)*d (A^fat)) {Lkys - Vk) ^ d{A''{y,),^j,{Lky. - Vk)t{Lky. - Vk)])^ 

2(Tk - s) ^ 2(Tk - s) >\ ^ ' 

where C is a positive constant. 

Proof. This one consists in using approximations y"^ solutions of ([5]) of process y solution of ([S]) . Thanks 
to Girsanov's theorem, we are able to obtain an equality for all bounded continuous function / 

E[f{x)G'{x)]=E[f{f)H^iyn] 

where /H^ is the density given by Girsanov's theorem. We want to prove that with a good choice for 
G"^ and 77^ , the lefthand member of the last inequality converges to the conditional expectation, and the 
righthand one converges to what appears in the Theorem [TJ 
We set for all z e R" 

fc=l ^ 

Then for all bounded continuous function / 

nfif)] = nf{x) exp{- KixtTdwt + ]^\\hl{xt)fdt}] 

We are looking for a different expression of the argument of the exponential function. We use Ito's 
formula iov Tk — eu < t < — e and use JT]) to get 

\\l3^{xt){LkXt - VkW\ _ 2{LkXt - Vk)* A'l{xt)Lkdxt , \\l3^{xt)(LkXt - VkW ^ 



Tk-t J Tk-t {Tk-tf Tk-t 

{LkXt - Vk)*d{A^{xt))iLkXt - Vk) d{Alj{x,), {LkX. - Vk)t{LkX, - Vk)j)^ 



Tk — t ^—^ Tk — t 
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The k term of (/if)*du't coming from that one in dxt is now isolated 

2{LkXt - Vk)*A^atdwt \\l3^{LkXt - '"kW _ _^ f \\l3tiLkXt ~ VkW \ _^ ruk 



Tk-t {Tk-ty V Tk~t J Tk-t 

2(LkXt~VkTA^btdt {LkXt~Vk)*dA'^{LkXt-Vk) '^{^Iv i^kX - VkULkX - Vk)j)f 

Tk-t Tk-t Tk-t 

(10) 

Since we have 

l!«tll a^ = Z^ fj^ _ l(T,-e,,T,-£)(i)dt 

and 

we obtain — 2(/if)*dwt — |pdt adding the terms given by ([TU|) . Finally, it leads us to a new expression 
for the density given by Girsanov's theorem 

^r,, m^_,{LkXn-,-Vk)r , \\P^r^_,^{LkXT,-e,-Vk)r 

E[/(y")] = E^exp I ^ ^ — + '-^ — 

k=l 

"^""^ {LkXt - Vk)*A%At nik 

Tk-t 2{Tk-t) 

{LkXt-Vk)*dA'l{LkXt-Vk) d(A,^^j-, {LkX - Vk)i{LkX - Vk)])^ -, 

.<km. 2(T,-t) > 

In an equivalent way, even if it means changing / 

nf{y'W]=nf{xW] (11) 

with 

„.,_.. . . \\Pk~eSLky?r,-e,-Vk)\\\ [^"-^ {Lkyl-VkYA%dt 



{y)= 11 Efe '7l(yT,-s)exp{2^ — — +/ 



(Lfcyf-^fc)*dA^(Lfcyf-^fc) ^ d(Afe^-, {Lky' - VkULky' - Vk)j)t ^ 
''^Tk-t) ^J-^^^^ 2(Tk-t) ^ 

and 

i/,^ = ||,y|(a;T,-£)exp{ 5^ — } (13) 



2£ 

fc=l 



where for all z £ M" 



7?|(z) = Jdet(A^^_,(z)) and C^ = n^" 



Now using it in the case where / = 1, we get formally 

E[f{if)^^] _ E[f{x)r 



the fact that this quantity is finite is given by Proposition [TJ The fact that the righthand term converges 
to the conditional expectation is given by Lemma IH] in the appendix. The proof relies essentially on the 
use of Aronson's estimates that provides gaussian bounds for transition probabilities. 
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The main difficulty of the proof consists in showing almost sure convergence and then uniform one 
for the Lp"^ . An obvious candidate for the limit is 



N 



= n '^fc ' '7fc(yTjexp{^ 



{yT^^eu){LkyTu-eu -Vk)\\ 



k=l 



2ek 



n 



{LkVt - Vk)*A^^Lkbt{yt)As 
Tk-t 



{LkVt - Vk)*d{A'^{yt)){Lkyt - Vk) 
2{Tk ~ t) 



{Lky. - Vk)t{Lky. - Vk)j)t -I 



l<i,j<.7nk 



2{Tk - t) 



Thanks to Lemma [TO] given in the appendix, Lp is well defined. As said before, we want to prove the 
following 

Lemma 3. There exists a decreasing sequence (ei)igN tending to such that 

lim E[|j(^'^' -(^||] =0 



Proof. The proof is decomposed into two main parts. First one aims at showing the almost sure conver- 
gence of (p^' . In second part we prove that E[(/3^*] tends to E[(/3]. Finally to conclude, we will use Scheffe's 
lemma. 

For almost sure convergence, we first use triangular inequality 

\v'{y') - ^{y)\ < l^'ivn - ^%y)\ + l^'iy) - 

The second one converges to 0, this is given by Lemma [TOl We now treat the term |(/3^(j/^) — ip'^{'y)\. 

N 



v'{y') 

ip^{y) 



11 «f^,^ A^'^'^X 



kJ{ VliVn-^e) ^- 2£fe 

r^"-' (Lky! - VkYAUyDLkhiy!) - {Lkyt - VkTAUytWtiyt) 



Tk-t 



dt 



{Lkyf - Vk)*A{A^i{yl)){Lkyl - Vk) - {Lkyt - Vk)*d{A'^{yt)){Lkyt - ffe) 



E 



2(Tfe - t) 

d(^t-(y")' ^^i^y' - ^k)7iLky' - Vk)j)t - d{Alj{y,), {Lky. - Vk)^{Lky. - Vk)j)^ 

2{Tk - t) 



We can write it respecting the order above 

^"{V"^) Notation 



N 



f^{y) 



k=\ 



According to Lemma[2]there exists a decreasing sequence (£i)igN tending to satisfying for all fc that 
y^^_^. converges almost surely to yxk- From this we obtain the fact that S^' converges almost surely to 
1 and T^' to using regularity of a. Then for all k 



m\ < 



Tk-e 
Tk-ek 



Lk{yt-ytYAl{yl)Lkht{yt) 



Tk-t 



{Lkyt - Vkr{AUy!)Lkbt{y!) - AUyt)Lkbt{yt)) 



Tk-t 



dt 



Since b and a are bounded we use Lemma [T] to get 



l*f.| < C 



^"-^ M-ytf 

Tk-Ck Tk - t 



dt 



Tk-e 



Tk-e 



(i + iogiog((rfc-0-i + , 

Tk-t 



-dt 



where C and C are positive random variables. Thanks to Lemma [21 up to an extracted subsequence 



lim 



\ Tk — t 



dt = Q 
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that leads us to convergence for all k of l^"^'] to 0. Now we use Identity ([M)) 
c't = — : Pt{yt)at + — - — - — qt{yt)dwt + —— —^rt{yt)dt 



T -t 



T -t 

\\Lyl~v\? 
T-t 



ft(yf)dt-^^f:^<Z*(2/f)d«;t 



where p, 5, and r are all Cl'^ functions. Hence using Lemmas [T] and [5] as above we obtain that 
lime._>o |0^' I = up to a subsequence. It remains to treat the term $|. Still using Identity (l34l) . 
Lemmas [T] and m we show that lim£._i.o 1$^'! — even if it means extracting once more a subsequence. 
We have obtained almost sure convergence of to Then we show the convergence of the expectations. 
For this we set a preliminary result 

Proposition 1. There exist two positive constants ci and C2 such that for all < e < Eq 

ci < C^EliP"] < C2 

Proof. We give an explicit expression 



/ g"(Ci, • ■ • , Cjv) n vUCk)~^ exp{ 
fc=i 



2e 



where is the density of {xti~e, ■ ■ ■ , ^Tn-e)- Under theorem's assumptions x is a strong Markov process, 
with positive transition density. For s,t € [0,T], we denote Ps,t{u,z) the density of x^'^ solution of ([J) 
initialized to be u at time s. Then thanks to Aronson's estimates there exist positive constant p,, A, M 
and A such that the density p satisfies for s < t 

_„ -A||z-u||^ _„ -A||z-u|p 

/i(t-s) 2e t-s < Ps,t{u,z) < M{t - s) 2e *-= 
Now using p we are able to write 

^''(Cl,- • -Xn) = Po.Tt-e{u,Cl) ■ ■ ■PTM-i-e.Tt.-eiCN-lXN) 

Then we apply Aronson's estimates and the fact that for all i,j the coordinate ^ is bounded by two 
positive constants. We obtain bounds for C^E[V'^] of the type 

^ ^11- - '|2 Ml/- „,||2 



A- 



-'C^J exp{^ 



-\\\LkCk-Vk\r AIICi-ull 



2e 



Ti -e 



A||Cfc - Cfc-l||% TT- , > 
fc=2 " *■ ^ fe=l 



where A is a positive constant large for the lower bound and small for the upper one. The integral can 
be interpreted as a gaussian expectation where 

/ Cr \ 

C2-C1 

\Cn - Cn~i) 

is a centered gaussian vector with covariance matrix 

/(Ti - £)/„ 



^ -2A 



V 











(T2-Ti)/„ 



(Tat - rjv_i)/„/ 



where /„ is the n-dimensional identity matrix. As a remark, in the sense of symmetric matrices, there 
exist two positive constants ci and C2 such that 



CiInu < < C2lNn 
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where I^n is the iVn-dimensional identity matrix. Thus the gaussian vector 

/Ci' 



admits for covariance matrix 
where 



G 



( In 

-In 


V 








■■. 

" Iji In / 



We still keep bounds for the covariance matrix 

CllNn < < C2lNn 

Now we can get bounds for C^E['(/i'^] with expectations of type 

N 



A^C^E[exp{-5:^^^^^p^}] 



fe=i 

where Xj- is a n-dimensional gaussian variable. Then we use Lemma [TT] given in the appendix to obtain 
the fact that 

M™i X • • • X M™" R 



ivi,-.-,VN) ^ C^nf=iA^E[exp{ 



2e 



'-}] 



is a gaussian density of a variable (£1X1 + y/jYi, . . . ,LnXn + y/jYN) where the Yfe are mfc-dimensional 
centered normalized gaussian vectors. Moreover the two families {Xk)k and (Yfc)fc are independant. 
Finally using bounds obtained above for we get the fact that for all < e < £0 



ci < C'E[i^'] < C2 



□ 



As a first consequence, thanks to identity (dH), E[(^'^] is finite so that g^sj is a density. We may also 
use Fatou's lemma to get 

E[ip] < liminf E[(p''] < C2 



e-i-O 



It takes more work to control limsup^_j.Q E[(^^]. 

Let J > be a large number, we introduce for all process {zt)te[o,T] and for all 1 < fc < the 
stopping time t| 



T^! = mi{tk <t<Tk-e : 
= mi{tk <t<Tk-e : 



1 



\\LkZt - VkW^ 
2(Tfc - 1) 



exp{-"^^;*-"^'''j} < J-i} 



2iTk - t) 



J 



} 



where _D is a positive constant such that Did < ^k- We know that such a constant exists according to 
assumptions on the function a. The tk are chosen to be real numbers contained in {Tk-i,Tk ~ e)- As a 
convention we set t| ~ Tk if the condition is empty. Let be the first of the t| such that the condition 
is non-empty 



inf{T| 



<Tk} 
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we set as convention r*^ = T if if for all k, — Tk- Even if it means changing / in Equation (ITT|) 
We recall that 

C'^=[[£ '^^P'f ^ J' 

k 

We now consider to be on set {r^ — '''k) 



Tk-i tk Tk — e Tk 

We decompose C^-iA^ into three parts as a product of three factors 

Fi = [[e 2 cxp{ } 

j<k 

-Uli r Wl^n^eiLkXn-e - Vk)\\\ 

F2=£ 2 exp{ i } 

Ze 

-Fs = 11 £ ^ exp{ ^ } 

j>k 

We are interessed in 

We now use Markov's property to get independance between Past and Future knowing Present (cf [3] 
see last chapter about conditional expectations) 

E[C^^^1,.=,=] = E[i^iE[F2^^3lr^=r||T|,a;..]] = E[FiE[F2l,.=,.|r|,x,.]E[f^3kr,d] (15) 

= nFinF2U^=rl\TlXr^]nFz\xT,-e]] (16) 

In order to study the factor E[F2li-==r| |=^r»] we introduce 

2(T. - ^) ^ 

For t e [Tk-i,Tk - e], we set zt = LuXt - Vk, pt = \\l3^{LkXt - Vk)\\. We recah that = a^LlA^ with 
Af = {I3^)*f3^ = (LkatLl)^^ . With respect to these notations, we have 

2 * Ak 

Pt = At Zt 

It is also easy to see that 

Azt = Lkbtdt + Lkatdwt 
and then d(z)t = LkatLldt = (A^)^^di. We use Ito's formula 

d(p?) = d{z;A'lzt) = 2ztA^dzt + zldA'lzt + ^ d{Al^,z^Zj)t + Wfedt 

Then 



Pt 



'2 2zfA*=dzt , , zldA^tZt , mfedt ^ Eij d(A,,j, z^Zj)* 



Tfe-i Tfe-< (Tfc-<)2 Tfe-t n-t Tk-t 
First using definitions of z, fi^ and A*^ we get 

zlA^tdzt = z^Afifedxt zlA'lLkatai^dxt 
= z*Mrai\dt + z*Mrdwt 
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This leads us to the existence of two bounded adapted processes r^^^ and r^^' defined on [Tfc_i,T/j — s] 
such that 

ztA'^dzt =ptrf^df+ptrf^d«;t 
In a same way we remark that there exist two bounded adapted processes r^^^ and r^^^ such that 

d{LkatLl 

we get even if it means changing r^^^ and r^'^^ 



d{LkatLl)=ri^'>dt + ri'''>dwt 



z^dA^Zt = ^r(d(Lfcati^)-i)^t =p?rf di + dtot 

Finally, we obtain existence of two bounded adapted processes r and r' such that 

2ptdwt , p'^dt pI dt ,p1+pt,, 

"^T-t = -T-T + + ^* T-^^"* +T-t+ '^^-T'^' 

Prom this we deduce that quadratic variation 

\Tk-t/ [Tk-ty 

Now we apply Ito's formula to the fimction 9 always for t € [Tfe_i, — e] 



2(rfe-t) 2 ' \Tk-tJ 8 ' \Tk-t 

We deduce from the three last equations after simplification of four terms that there exists a martingale 
M and a bounded adapted process r" both defined on [Tfc_i,Tfe — s] such that 



For any 1] > 0, functions x t-^ e ''2|^|'"form=l,2,3,4 are all bounded, then there exists a constant 

,4 I ^3 



Cn such that 



This gives us the existence of a bounded adapted process tt defined on [Tk-i,Tk — e] that allows us to 
write 

d0t = dMt + e]-\Tk - ty'^atdt 

with h = We now integrate it for t € {T^,Tk — e] 

0t =- + Mt - Mrt + / 7r,0i-''(Tfe - s)-^ds 

This leads to the following 

netiri<t] <J-' + 7t[ E[0,i,.<,]i-''(Tfe - sy'^ds 

Jtk 

where tt = sup^ \ns\. So E[61x|<t] is bounded by the solution u of 

dus = au\-\Tk - s)-^ds, ut^ = 
and this equation has an explicit solution 



Ut 



1_ 
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where Ck is a positive constant. Then for all t G [Tk^i — Tk — s\ 
In particular when t ~ Tk — e 

We now come back to equation (IT5|) . we get a first bound 

E[C^V'lr.=r,d < {ck{Tk - tk)'-'' + J-TE[i^iE[F3|a;T,-e]] 
In order to treat the factor E[F3|a;Tfc-£] we use Aronson's estimate to get 

E[F3\^T,-e] < G / _[]_ — exp{ i }- -—expi-AjLl I^jdQ 

where G is a positive constant. We just have to use Lemma [TT] given in the appendix to obtain an 
positive constant upper bound. The same Lemma [TT] brings us a positive constant upper bound for 
E[Fi]. Finally the inequation we get from equation (TTSl) is the following 

ncrU'^rf] < G{ck{Tk - tk)'-'' + J-"}" 
where G is a positive constant. From this we deduce 

E[G^V"lT>r^] =5]E[G^V'lre=r|] < G max { (cfc (Tfc - ifc ) + J-")" } 

k 

According to this last result and using the lower bound of G'^E[';/'^] given by Proposition[T]we finally have 

where G is positive constant. So using inequality (jlip we obtain 

^1^^ > 1 - Gmax{(c.(r. - tk)'-" + J"")"} (17) 

Moreover the family (ip'^lT=T')e is uniformly integrable. Indeed by definition of we can get upper 
bounds depending on J for the different factors in Expression ([T^ of (p^ or (jTT]) of (p, for all < e < 1 



e T\ e( e ^ r f^"'" {LkVl - Vk)* Lkbtdt {LkVl - Vk)* dA^ijLkVl - Vk) 

V- d(A^_^-, [Lkv' - Vk);{LkV' - Vk)])t ||/3T.-.,(^/T,-sJ(-^fc2/T.-5. " 
^ 2(Tfc - 1) 2ek 

L fact 'V*' — </'"• We recall that b and ct are bounded so is rj. Then on {r^ = T} 



Tk-t 



C 

< 



log 



J 



which is an integrable quantity in Tk, and G is a positive constant depending on the choice of b and a. 
A same method gives an upper bound for the terms where quadratic variation appears. For the terms 
in dAf , we decompose with respect to integrals with respect to dt and dwt 

iLkvi-vknAHLkyi-vk) ^ MM_L^,.(,e),^^ ii^f 

^Ufe ~ t) J-k ~t J-k — t 

where and are bounded adapted functions. Then for fixed J, there exists a constant K such that 

^n... < Gnexp{ r - i "^,y-;;f ik^i^dn 

fc M-efc ^fe-r ^ {-ik-ty 

where G is a positive constant. Let us recall the following lemma (cf J) p. 198) 
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Lemma 4 (Novikov). Let (Mt)tgR be a continuous local martingale, we set for all t 



If 



then we have 



Zt = exp{Mt-i(M)t} 



E[exp{i(Af)t}] < +00 



Let us remark that for all p > 



ex 



log 



\ H 2 

J 



dt} < CP 



where C is a positive constant. Thus, we apply Novikov's lemma to get uniform integrability. Then we 
take the liminf£_j.o and use Lebesgue's theorem to obtain 

^^'^^^=-'^ > 1 iVmax {(cfe(rfe - e t,f-'^ + J-')}^ (18) 



lim supg_j.Q E[((5'^] k 

Now lT=r= converges almost surely to 1 as J tends to infinity. We are able to say after making the 
tend to Tfe that 

limsupE[<y9''] < E[(^] (19) 

e->-0 

We finish the proof by Scheffe's Lemma (cf [3] p. 36) 

Lemma 5 (Scheffe). Let {fn)n<£N be a sequence of positive functions converging to f , moreover we 
suppose that 

lim E[/„] = E[/] < 00 
then the sequence (/„) converges to f inlJ-. 

□ 

Finally we conclude thanks to Lemmas [9] and [3l □ 
Case where b is unbounded 

Suppose now that b is locally Lipschitz with respect to x and is locally bounded. Moreover the SDE ([1]) 
admits a strong solution. We use a Girsanov theorem to reduce the problem to the case of a bounded 
drift. 

We recall the Girsanov theorem for unbounded drifts introduced in |3] 

Theorem 2. Let b, h and a be measurable functions from x R" to M", and R"^'' locally Lipschitz 
with respect to x; consider the following SDE's 

dxt = bt{xt)dt + at{xt)dwt, 

<iyt = {h{yt) + (Jt{yt)ht{yt))dt + crt(yt)dwt, 

xo = yo 

on the finite interval [0,r]. We assume the existence of strong solution for each equation. We assume 
in addition that h is bounded on compact sets. Then the Girsanov formula holds: for any non negative 
B Orel function f defined on C([0, T], M"), one has 

E[f{y,w'')]=E[f{x,w)exp{ f hl{xt)dwt-\[ \\ht{xtWdt}] 

Jo ^ Jo 

T 1 '■^ 



E[fix,w)]=E[fiy,w'')exp{- I hl{vt)dwt - ^- j \\ht{ytWdt}] 
where — wt + /q hs{ys)ds. 
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Thanks to both Theorems [T] and [5] we obtain 

Theorem 3. Suppose a and a"^ to be C^'^ -functions. Assume that b is a locally Lipschitz with respect 
to X and locally bounded function. Let y be the solution of 

dyt = k{yt)dt + at{yt)m -J2Myt)f3Hyt) 'J!'_J' lin-e,,n){t)dt 

k—1 

where b satisfies the assumptions of Theorem^ 
Then for any bounded continuous function f 

E[f{x)\{LkXTk ^Vk)l<k<N] 



CI 



^ II Rk (T.,.^ M|2 



||/3T,-efc(^fcyT,-e. - Vfc)lr , [ (Lkys - Vk)* Lkb,{ys)ds 



{Lkys - vk)*d (A'^iyt)) (Lkys ~ Vk) ^ d{A''{y,),j, {Luy. - Vk)i{Lky. - Vk)j)^ 



2(Tk - s) ^ 2{Tk - s) 

^ 1, 



+ / bl{yt)at{yt)-^Ay,-^\at{yt)-^b,{y,)fdt) 





where C is a positive constant and b — b — b. 
Proof. Let x be the solution of 

dxt ^ bt{xt)dt + at{xt)dwt, xq ^ u 
Then from Theorem [21 for any bounded continuous function / and g 

E[f{x)g{L,XT„. . . , Lnxt,)] = E[f{xt)g{L,XT„ . . .,LNXT„)e^o 6:(£.)a*(«0-Mx.-i||a.(s.)-i6.toOll^dt] 

k 

It remains to apply Theorem [TJ □ 

Appendix 

Lemma 6. Let us consider Equation l^jj with random initial condition u on [0,T] with N = 1 which 
means only one observation time in T . 

dyt = bt{yt)dt + at{yt)wt - ft(2/t) ^y l(T-ei,T)dt, yo = u 

Then this equation admits a unique solution on [0,T). Moreover \\L(yt~ui)\\'^ < C{uj){T ~t) loglog[(T — 
t)^^ + e] a.s., where C is a positive random variable. 

Proof. We recall that parameters b and a are locally Lipschitz functions. So that the equation admits a 
unique solution on both intervals [0, T — £i\ and (T — ei, T) and so on [0, T). Moreover thanks to Ito's 
formula, on (T — Si, T) 

d^fc^ = (T - t)-'L[hdt + a^dm - Pt'-^dt] + L^l^dt 
T-t ^ ' ^ T~t ^ (T-ty 

then using we have LPt = L so that 
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For all 1 < i < 71 the process {{Jq{T — s) ^(Tsiys)dws)i,t > 0} is a continuous local martingale whose 
quadratic variation Tt = J* X]j=i(^ ~ s)~-^as{ys)i.jds satisfies linit_j.T Tt = +oo and Tt < jf-^ where c is 
a positive constant. Hence we just have to apply the Dambis-Dubins-Schwarz theorem that gives us the 
existence of a Brownian motion such that 



I {T - s) ^as{ys)dws I 
/ , 



The law of iterated logarithm allows us to conclude. 



□ 



Lemma 7. Let us consider Equation with random initial condition u on [0,T] with N = 1 which 
means only one observation time in T 



Then for all s < t <T , 



and 



dyt = bt{yt)dt - Ptiyt)^Y^j-dt, yo = u 



T-t 



< c(l + VT^E[\\L(u - ui]\\^]) 



n\ys - ytf] < c{t - s){i + Vt~:s¥.[\\l{u - ui)f]) 

where c and C are positive constants depending on T , ei, bounds for b and a. 
Proof Thanks to Identity Q, on (T - £i, T) 

dL{yt — U\) ~ LSbidt + atdwi\ — —dt 



T-t 



Thus 



d{\\Liyt - = 2iyt - m)* L* L[btdt + atdwt] - 2^^^^^^^^dt + TiiLatL*)dt 



T-t 



where the function Tr gives the sum of all diagonal terms. Finally 



\\Liyt - ui 



Setting Et=E 



T-t 



\L{yt-ui 



T-t 



= 2 (^i^L*L[5,dt + a,dm] + '^^^^^dt 

T-t L t T t tj 

, since b and cr are bounded, we get 



\Liyt-u^W 
[T-ty 



dt 



E't<Ci 



/Et + 1 
T-t 



T-t 



where Ci is a positive constant depending on ||6||oo and \\a 

E[<{T-tr^ 



= {T-t)-\C-^) 



where C = Ci + ^. Thus 



Et-2C\' 



thanks to 



Hence 



Et-2C 



J 2(T-t)! 



Et — 2C ^ Ex-ei ~ ^C" 

y/T^t ~ 



< 



that can be written 



Et<2C 



T-t 



{3 



■T-El 



2C) 



(20) 
(21) 



(22) 



(23) 
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Similarly for t < T — ei, Inequality [21] becomes 

Et<C\Eo + l)exp{ — } 

El 

where C is a positive constant only depending on T and bounds of b and a. That gives us ([101) • 
By definition for s,t G {T — ei,T) we have 



Vs - yt 



b-rdr + Grdw-r — Pr^^ -dr 



T-T 

Since b and a are bounded functions, using Minkovski's inequality 

nWVs -ytf]'^ < n\\ J' Mrf ]^ + E[|| J' ardwrW']^ + E[|| j'^ P.^^drf]^ 
Thanks to Doob's inequality (see e.g. [5] p. 170) we get 



(24) 



E[|| ^ brdrf] + E[|| ^ fx^duviri < C2(i - s) 



where C2 — Cf is the square of the constant introduced above. In order to treat the last term in ([M]). 
we beforehand give a property 

Proposition 2. Let f be a real-valued process defined on a segment [a,b], then 



E[{ [' fsds)']i < fE[f^]ids 

J a J a 



Proof. Indeed 



so that 



E[(/ fsdsf]<E[{f \f,\dsf] = 2E[f f \f^\du\fs\ds]=2 f f E[\fufs\]duds 
E[( ( fsdsf\ < 2 r rE[(/„)2]^E[(/,)2]^duds = ( [\[f^]^dsf 

J a J a J a J a 



□ 



Thanks to Proposition [H assumptions Q on matrix P and result (j20p 

\ 2 



Pr 



T-T 



dr 



< 



\\\\PAyr-u,)r 



(T-T) 



ds 



< c(l + Vt^E[||L(u-ui)||2])( 



dr 



-.y ^4ciT-s){l + VT^E[\\L{u-ui)\\^]) 



Finally using {a + b + cy < 3{a^ + b^ + c^), we have on (T - ei,T) 

E[\\y,~ytf] < {C2AAc){T~s){2 + Vt~sE[\\L{u-ui)\\^]) 
Using Doob's inequality for s, i e [0, T — ei] with s < t 



E[\\y,-ytf]<C2{t-s) 



this gives the second result ([21 



□ 



Lemma 8. Let y and z be two bridges, solutions of (3^ with N — 1 and different initializations. Then, 
there exist two constants C > and < a < 1 such that for all t € [0, T] 



\yt^ztf]<CE[\\ya-z^fr 



(25) 
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Proof. Using their definition, we have on [T — Si, T] 

-/ N f \ u f \^A-^ , \ t \ f \^A~ Ptiyt){yt - Ml) - Pt{zt){zt - Ul) 
d(yt - zt) = [ot{yt) - ot{ztmt + [at{yt) - (Tt{zt)\Awt Y~t 

Thus 

d||yt - ztll^ = 2{yt - zt)* [[bt{vt) - bt{zt)]dt + [at{yt) - (Tt(zt)]dwt 



Pt{yt){yt - Ul) - Pt{zt) {zt - Ul) ^ 2 /ofi^ 
dtj + 2^(crt(2/t) - at[zt)),^jdt (26) 



In a same way, on [0, T — Si] we obtain 

d\\yt - ztf = 2{yt - Zt)* [[bt{yt) - bt{zt)\dt + [at{yt) - cjt{zt)]dwt] + Y.{<Jt{yt) - <^t{zt))ljdt (27) 

We denote Et = ¥\\\yt - zt|p]. We decompose the interval [0,r] into [0,r - h] and (T - h,T] with a 
parameter h that will be chosen later. On [0, T — h] with respect to both precedent Inequations ([25]) and 
((27)) we get by using regularity of b and a 

where Ci and C2 arc positive constants depending on T, b and cr. We use Gronwall's lemma to obtain 
For the other part (T - /i, T], we use (US]), (H?]) and ^ to get 



where C3 and C4 are positive constants depending on T, b and cr. Then 

\og{Et + 1) - iog(i;T-,» + 1) < C4(V^ - Vr^) < CiVh 

Finally, on [0,r] 

Et < Eo (^^^ " (1 + e"^''^) + e^=^ - I < Eo (J^^ \l + K) + C^kVH 

— 2 _J — zu 

which minimizes the last member above. Hence 



Et < CE^''^+' 

where C is a positive constant depending on b, a, P and T. □ 

Proof of Lemma\^ . We now consider an interval of type [Tfc_i, Tk). We introduce a process y^ solution 
on this interval for the Bridge Equation ([3]) initialized at time Tk-i by the value yf^^ ^. A picture to 
visualize what is going on is given by Figure [T] in page 1171 We use this new process y'"' to write 

EM ~ ytf] < 2EM - in + 2E[|!y,'= - ytf] (28) 

We will study both terms separately. 

For the first one, on [Tk-i,Tk — e) the term \\y^ — y^\\ is a.s. and on [Tk — £, Tk), we can reduce the 
study to that of ||a;t — ytW^ with a same initialization yj'^_^ at time T^ — s. 

E[E[\\yi ' ytr\yk-e]] < '^Knhk^e - ytr\yk-e]] + 2E[E[\\yi - vk-A'lyk-e]] 
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Figure 1: Illustration of the three different dynamics considered 




Let us recall 

dx = bt{xt)dt + at{xt)dwt ^ 
dy = btiyt)dt - Ptiyt)^^^;^ + atiyt)dwt © 

(T) First, the three processes follow the dynamics of the initial diffusion ^ with a different initialization 
for y. 

(D Now, the three processes follow the dynamics of the bridge (jS]), that means that the correction term 
operates and forces these processes to get closer to the observation 

(D At the end, only the processes y and y'^ go on following the dynamics of the bridge Q and both 
tends to the obervation, while follows the initial dynamics ^ 
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We then use Lemma (O given in the appendix and classical technique (see e.g. JB^ p. 170) to obtain 
upper bounds 

E[EM - ytf\yk~e] < + V~eE[\\Lkiy'T,_^ - u,)f]) (29) 
where c is a positive constant. Now in order to treat the remaining term we use Lemma [8] 

Finally, on [Tk-i,Tk) 

nU - ytf] < c' [e{l + V^E[\\L,{y^^^_^ - u,)^]) + E[\\f^^_^ - yT,_,\\T] (30) 

where < a < 1 and c' is a positive constant only depending on T, bounds b and a. We show by 
induction that there exists some constant C such that for all 1 < k < N 

E[bf^,-2/TjP]<C'=£"'="' (31) 

The base case is given by Equation (I^Hl) . Indeed on [0,Ti] processes y^ and y are indistinguishable since 
they have a same initialization at time 0. Suppose now for some k that Inequality (|3ip holds. We now 
use Equation ((50)1 to get 

n\\yk+^ - 2/T,+JP] < c' [e(l + yTeEiWL.+.if^^ - u^+Of]) + E[||y|,^ ~ yT.f]"] 

Let us recall that LkyXk = -^fe^fc hence 

E[\\Lk+iiy'T, - Uk+iW] < E[\\Lk+i{y'T, - yn)f] + \\Lk+i{uk - Uk+iW < Ck{l + E[||2/f^, - yn\\^]) 

where Ck is a positive constant depending on Lfc+i, Uk and Uk+i- That gives us thanks to the induction 
hypothesis 

IE[||yT.+, - < 4£(1 + V^IE[||yf.^ - 2/TjP]) + IE[||y!., - yrJP] 

where C is a positive constant. This concludes the proof. □ 
Lemma 9. Let (tk.q) i<k<N be a sequence such that tk^q G {Tk~i,Tk) and for all k, (tk^q)q is an 

l<q<Mk 

increasing sequence. Then for all bounded continuous function g 

^pj^ ^ E[g{xt,_,, . . . ,Xt„_^,J\(LkXTk =Vk)i<k<N] 



Proof. Let us recall 

'2s 



C = lie 2 77|(xT,-s)exp{ } 



k=l 

where for all z £ 



r,|(z) = ^det(A^^_,(z)) 



Let introduce Aronson's estimates (see e.g. [I], [8] or [2]) that gives bounds for the transition density. If 
Ps,t{u, .) (with s < t) is the density of Xt knowing that Xg = u, we have for all z 

//(t-s) 2e < Ps,t(M, z) < Af (i - s) 2e *-= 

The transition densities allow to expand the density q"^ of [xt-^ n ■ ■ ■ , Xt,^ i xti-6i ■ • • i sjtjv-e) 

9" (2:1,1 ' • ■ • ' 2:tjv^„^ ,Cl, ■ • ■ , Cat) = PO,ti,i(u, Zi^i) . . .pt^ j,,^^,Ti-e(^0,Mo, Cl) ■ • ■ PtN ^Mj^ ,Tn -eizN,MN Xn) 
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Then we set for e > 

= / 9(zi,i,---,zt,^,M„W{zi,i,---,Zt„,M^Xi,---,CN)Y[dzj 
This application is continuous according to Aronson's estimates. From this expression it comes 



We recah that the rows of each matrix Lk form an orthonormal family. We now complete arbitrarily 
each family into an orthonormal basis of M". We denote Pk an arbitrary matrix whose first rows are 
given by ifc. Then we make a basis change with respect to those matrices Pk for each Cfc. Thus 



^l{P^\^. . . . , P],\m) n vl{Pk'Ck) eM- '^^X ^}da 

k 

denoting C^-' the vector composed by the coordinates from i"^ to j**^ one of C,k- We now make a second 
change 

r ek"'''^^eek"'''+vk 

i — ^k 

So that 

We now use Aronson's estimates and Lemma [TT] to get an integrable uniform upper bound for q"^ when 
< e < Eq. Thanks to Lebesgue's theorem we obtain the convergence for the last term 

/• II qk f l:mfc ||2 

n / (p.-' (c--) >->(^r^ (c--)) n^^-p{-^^^K. 



We then integrate with respect to the £,1'™''' 



n - J {p.-' (c--) ' • ■ • ' (A-^ {cH) n^c-^^^ 

k 

Finally 

lim 4 = " ^ ) ' ( _V2~ " (32) 



/<i>°l 












[Pi' 1 




\,---APi'\ 




)n.der+i^" 



We conclude thanks to the Bayes formula. □ 

Lemma 10. Let y be solution of Equation Then almost surely for all 1 < k < N the following 

integral are absolutely convergent 

{LkVt - Vk)* A^{yt)bt{yt)dt ^ [Lkyt - Vk)*dA'^{yt){Lkyt - Vk) 

n-ek Tk-t ^Tk- t) 

yr^ d{Al^{y.), {Lky. - Vk)t{Lky. - Vk)])^ 
^ 2{Tk-t) ^-^-^^ 
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Proof. We reduce the study without loss of generahty to that of 

dyt = bt{yt)dt + at{yt)dwt - cTt{yt)fit{vt) ^^ _ ^ l{T-ei,T){t)dt 

We then treat integrabiUty for each term. 

For the first term, since b and {3 are bounded, we use Lemma [T] to get 



Lyt - V 



T-t 



<C\ 



^ loglog((r-t)-i+e) 
T-t 



where C is a positive random variable. Now for all positive a, we have log log a; < x". Then for a small 
enough, we obtain integrability of righthandside. 

For the second term in we recall that for all z we have At{z) — (3t{z)* (3t{z) — {Lat{z)L*)^^ , 
hence 

dAt = ptdt + qtdwt + rt _ ^ dt 
where p, q and r are bounded adapted processes. So that even if it means changing p q and r 

{Lyt-vydAt{Lyt-v) \\Lyt-vf , \\Lyt-v\? . , WLyt-vf 

— -ptdt-\ — — - — qtdwt-\ — n dt 



T-t 



T-t ^ T-t 

'^yt-v\\^ \\Lyt-v\\'^ 

T-t ' {T-ty^ 



{T-tf 



(34) 



Using Lemma [TJ we obtain that the quantities , ^^"^^'_^"| and are integrable in a left 
neighboorhood of T. 

For the last term in (j33p . we use Ito's formula and the fact that Lat{z)l3t{z) — Id, so that on 
iT-ek,Tk) 



d{Lyt -v) = L[btdt + (itdwt] - ^^j^—^dt 



Hence 



d{A^^J, {Ly, - v),{Ly, - v)j)^ < \\Lyt - v\\ptdt 
where p is the same bounded adapted process given above. Finally 

d{Aij,{Ly, - v)i{Ly, - v)j)^ ^ \\Lyt - ^WPt 



T-t 



T-t 



even if it means changing p, and this last term is integrable. 



□ 



Lemma 11. Let {Zj)Kj<,K be a family of random ruj -dimensional variables and let [g.j 
^)i<j<K be a family of densities. Then the function 



^ E 



is the density of the family [Vj = Wj + Zj)i<cj<K where each of the Wj whose law is given by gj is 
independent with respect to the {Zj)i<cj<K o,nd (Wk)k^j. 

Proof. Let / be a bounded continuous function 



Y[gj{Zj-Vj) 



Y[d^ 



Then we make the change of variables Wj = Zj — Vj for all j 



l[dv,=E 



f{iwj + Zj)j)Y[9jiwj] 



Y[d^ 



Y[dwj 



Hence 



E 



f{{wj + Z^)^)'[{gj{wj) 



l[dw,=E[f{{W,+Z,),)] 



where Wj admits gj as density. 



□ 
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